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Abstract 

The role of electronic interactions in the level structure of semiconductor 
quantum dots is analyzed in terms of the correspondence to the integrability 
of a classical system that models these structures. We find that an other- 
wise simple system is made strongly non-integrable in the classical regime 
by the introduction of particle interactions. In particular we present a two- 
particle classical system contained in a d-dimensional billiard with hard walls. 
Similarly, a corresponding two-dimensional quantum dot problem with three 
particles is shown to have interesting spectral properties as function of the 
interaction strength and applied magnetic fields. 
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I. INTRODUCTION 



The role of electronic interactions on the observed properties of semiconductor quan- 
tum dots has increasingly been found to be of vital importance, as many of the papers 
in these proceedings and elsewhere illustrate. More and more, transport and capacitance 
experiments, as well as detailed studies of far-infrared response, require the consideration 
of interactions in order to understand the observed experimental features. Since a number 
of experiments explore specifically the electronic level spectrum, it is of interest to contrast 
these studies with a detailed theoretical analysis of the role of interactions. We present here 
an analysis of the energy level statistics of a quantum dot system as function of variable 
interaction strength, and then as function of an applied magnetic field. 

As the level spectrum is considered in detail, the correspondence with the dynamical 
integrability of a classical system is also of interest. It is anticipated that as the Coulomb 
interaction is introduced, the dynamics would in general become chaotic and this would in 
turn be reflected in various statistical properties of the level spectrum. It has been known 
for some time now that as a classical system becomes non-integrable [[|], the corresponding 
quantum system exhibits a level spacing distribution given by the 'Wigner surmise' derived 
in the context of random matrix theory (RMT) 0. In fact, this behavior has been verified 
in a number of theoretical and experimental systems, although typically the classical non- 
integrability is due to disorder or boundary (geometrical) effects. Examples discussed in 
the literature include small disordered metallic particles |§ , particles moving in a variety of 
'stadia' and in two-dimensional antidot arrays RMT also has been used to make 

definite predictions on the statistical distribution of Coulomb blockade peak amplitudes . 
This behavior has in fact been shown recently to be a good description of quantum dot 
systems in beautiful experimental realizations fL0 |. One should also mention that recent 
work on interacting systems, some without an obvious classical counterpart, has also shown 
that these exhibit the level repulsion characteristic of quantum versions of non-integrable 



classical systems [pTj. Moreover, recent analysis of the level spectrum of excitons in quantum 
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wells (via photoluminescence excitation spectroscopy) has shown evidence of level repulsion 



g. 

A level structure described by RMT has been sought recently in the charging (or addition) 



spectrum of quantum dots ||13|| . Detailed analysis of experimental charging energies, after 
proper subtraction, would be expected to yield the single-particle-like excitation spectrum 
described by RMT. Unfortunately, it appears that the extraction of this excitation spectrum 
is obscured by the systematic shift in the charging energies, and the resulting level spacing 
is described by a normal distribution, rather than by the RMT functions. This would then 
suggest that a more direct measurement of the excitation spectrum (via perhaps nonlinear 
tunneling currents) would be desirable (although clearly difficult experimentally beyond 



the first few excitations |1J])- One should also mention here that it is believed that the 
nonlinear transport experiments explore mainly the excitations of the center of mass of 
the system (in the typically parabolic quantum dots) |L5|], due to the strong electronic 
correlations suppressing most of the tunneling 'channels'. This prevalence of the center of 
mass excitations is however expected theoretically to diminish as the energy of the excitation 
increases This regime would be reached only as the bias voltage is raised in transport 
experiments, and makes it then difficult to achieve in practice. We hope, however, that the 
results presented here would motivate more experimental work in this direction. 



II. A CLASSICAL DOT 

We have studied the effects of particle interaction on the classical integrability of a system 
of two masses moving inside a d- dimensional billiard. We find, in general, that the motion 
is strongly chaotic (typically exhibiting 'soft' chaos, with a mixed dynamics, where regions 
of phase space are still periodic or quasi-periodic), and with a strong dependence on the 
characteristic interaction length and strength. As perhaps one of the simplest examples (see 
Ref. |16|] for a description of the general case), consider two particles of equal masses moving 
in a ID box defined in the interval q : (— |, |) (we measure all lengths in terms of the box 
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size). The particles are assumed to interact via a screened potential V(qi, 52) = exp(— A|gi — 
Q2\)/\qi ~ fel j where A is the inverse screening length. Notice that this potential goes to a 
hard-core S— function when A —>■ 00, and the particles behave then as non- interacting but 
impenetrable points. In that case, the dynamics can be integrable in special cases of the 
mass ratios [T7]|. In this sense, A plays the role of a perturbation parameter which changes 



the degree of integrability of the system, as it determines the effective 'radius' of the particles 
for a given total energy. Direct calculations show that this is indeed the case. 

The problem could in principle be solved by direct integration of the equations of motion, 
but we find convenient to transform it to a set of center-of-mass and relative coordinates 
R — (li + ^2) / M , and r = c^ — <?i, respectively, where the total mass M = 2, and the reduced 
mass fi = 1/2. These equations define a two-dimensional space of coordinates p = (r, R). In 
this space, we have a new set of equations for the boundary of the billiard, say Fj(p). The 
Hamilton equations are transformed to r = p/p, R = P/M, and 

p = -V r V(r) + £ Mp, P) 5[F 3 (r, R)} 
3 

P = Y,B J (p,P)m(r,R)}- (1) 
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The functions Aj and Bj describe the change in the momenta p and P, due to the bounce 
on the j-th wall. Notice that these equations describe the motion of a single 'hyperparticle' 
in the two-dimensional p-space. This hyperbilliard description can be generalized to any 



number of dimensions 16 



Notice that bounces of the hyperparticle in the hyperbilliard correspond to bounces of the 
masses in the real/original dot. The walls of the billiard cause the breaking of translational 
symmetry of the system, and as a consequence, the center-of-mass (CM) momentum is no 
longer a constant of motion. In the case of non-interacting and equal-mass particles, the 
changes in the CM momentum are determined only by the geometry of the billiard. In 
our case, however, the interaction couples the CM and relative momenta after each bounce, 
which in turn depend on the momenta of each of the original masses. 

Dynamical map. — We should also mention that these equations in p-space provide an 



interesting description which is also extremely useful: In between bounces, the hyper-particle 
moves freely along the i?-axis, whereas the interaction acts only along r. These two motions 
are independent, and only become correlated at each bounce, when the different momentum 
components are changed, while keeping the total energy constant. Understanding this fact 
allows one to describe the motion in terms of a dynamical map connecting the different 
bounces. If the coordinate of the hyperparticle is p n = (r n , R n ) at the time of the n-th 
bounce, then the time spent until the next bounce on the j-th wall is obviously the same 
along both components, and one can then write r r (p n , p n +i) = tr(Pu, Pn+i)- Since the CM 
motion is that of a free particle, except for the collisions with the walls, tr can be calculated 
simply. Similarly, for a pure Coulomb potential (A = 0), r r can be calculated analytically, 
so that the previous equation can be written in a more explicit form, r{T e (p n ), T e (p„ +1 )} =| 
(R n +i — R n )/(P/M) |, where r is now the time elapsed going from p n to p n+ i, expressed 
in terms of the time T e (p) spent by the particle from the turning point to p, T t (p) = 
if + (2^) ^ costT^re) 1 / 2 . Here e = E — P 2 /2M, represents the energy left for the relative 
motion, as E is the total energy of the two-mass system. For a weakly screened potential, 
A < 1, we can expand V(r) to first order and obtain a similar expression, where E is scaled 
to E — A. What follows, after these definitions is to characterize all the different possible 
trajectories in the triangular region in p-space. A simple algorithm can then be obtained to 
determine the Poincare surfaces of section. This nontrivial (and clearly nonlinear) algebraic 
map provides then a full description of the dynamics. Its use (in lieu of the direct integration 
of the equations of motion) simplifies calculations a great deal, and allows one to better 
characterize the system, as we describe below. 

Poincare sections. — In order to characterize the motion of the two interacting parti- 
cles (or the hyperparticle with its many degrees of freedom), we explore Poincare sections of 
the resulting phase space. Having a four-dimensional space in this case (r, R,p, P), we select 
to show those where one of the particles is at one end of the box (notice that the energy 
is a constant of motion here). Figure 1 shows a typical Poincare section for A = 0.1 and 
E — 1.5, obtained using the map described above (and which is virtually identical to the 
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one obtained directly from integration of the equations of motion ||16|| ). Notice that chaotic 
trajectories nearly fill the available phase space (for this given total energy). Increasing 
values of A give chaotic orbits that fill more of the available phase space. Moreover, there 
are also a number of islands of stability, as expected from the KAM theorem [|1], near the 
fixed point corresponding to periodic symmetric motion in the non-interacting system (and 
indicated in the figure with a cross on the right axis). In fact, for all values of E and A, the 
initial condition where qi(0) = —92(0) = |, and pi(0) = —p2(0) = \J E — V(l) gives rise to a 
periodic orbit. Other islands also appear purely due to the interaction, as can be seen near 
q2 = 0.2, P2 = 0, and are associated with a sort of correlated motion of the two masses. 

General values of E and A give rise to this type of mixed dynamics, with chaotic and 
regular trajectories sharing the available phase space. This increasing degree of 'soft chaos' 
in the system (as A increases, for example), should be reflected in the level statistics of 
the corresponding quantum mechanical system, as we describe explicitly elsewhere in these 
proceedings. One naturally expects that this effect of interactions — turning a regular 
system into a non-integrable one — would be rather pervasive, regardless of the type of 
particle confinement and details of the interactions. In the following section we illustrate 
this effect in a somewhat different quantum mechanical model of a quantum dot. 



III. A QUANTUM DOT 

Here, we will use a model of a quantum dot which has been very successful in the 
description of recent experiments in these structures f6,15|,18l. The quantum dot is modeled 



as a parabolic potential well with circular symmetry, and the frequency (or curvature) is 
chosen to fit characteristic single-particle excitation energies in these devices. In this case, 
the few-particle problem, which includes fully the effect of interactions, can be solved quite 
accurately (numerically 'exactly') even in the presence of magnetic fields. The approach 
is based on a canonical (Jacobi) transformation to a set of auxiliary harmonic oscillator 
generating operators which allow one to write the interaction matrix elements in a closed 
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analytical form, easily calculable ||19|| . This, together with the ability to separate the center 
of mass (CM) degree of freedom from the 'relative' ones (a feature almost exclusive to the 
parabolic confinement potential), allows one to completely characterize and solve for the 
spectrum of up to three electrons in this parabolic well. Correspondingly, the spectrum can 
be characterized by E = e(N C m, M C m) + -^rei, where the first term gives the CM manifold, 
and E ve \ is obtained from the relative-motion part of the Hamiltonian. In the case of three 
particles, this latter part has diagonal elements given by e(iV_,M_) + e(N + , M + ), which 
are mixed by the Coulomb interaction (taken here as given by Vfaj) = e 2 /(erij), with e 
a background dielectric constant which would be assumed variable in the results section 
below). Notice that in all these expressions, 

s(N, M) = hQ + + m - (M + -) , (2) 

arise from the auxiliary harmonic oscillators introduced in the canonical transformation 
I9[, with fi± = \/4fi;j + &a ^ w c /2- Here, characterizes the single-particle harmonic 



confinement potential of the quantum dot (with a typical value of HQq ~ 1~2 meV), and 
u c = eB/mc is the cyclotron frequency of the electron. Due to the rotational symmetry 
of the system, the Coulomb interaction only couples states with the same relative angular 
momentum, given by L Te \ = h(M + + M_ — N + — AL). This is important, for apart from 
making the calculation simpler, it also allows us to carry out matrix diagonalizations with 
extremely high accuracy and for very many eigenstates (typically one to two thousand levels). 
This complete convergence is obviously important if one is interested in the analysis of the 
level spacings, as we now proceed to do. 

Level statistics results. — The statistical analysis of the levels obtained as described 
above for three particles can be carried out according to the typical prescriptions in the 
literature. The level spectrum is used to extract a slowly-varying density of states with a 
smooth energy dependence, characteristic of the system at hand. This unfolding procedure 
then leaves one to study the structure of the level fluctuations of the spectrum, on which a 
number of statistical tests and statements can be made The unfolding here is performed 



by defining the 'staircase' function N(E) which gives the cumulative number of states below 
E, and then fitting this to a smooth polynomial (typically of fourth degree), N(E). The 
'linearized' or unfolded level sequence is then obtained from X{ = N(Ei), where Ei is the 
original sequence. This process, although not unique, gives similar results to other unfolding 
procedures (see jTI]] for a good discussion). 

One can perform a number of statistical analyses. Here we focus on the nearest-neighbor 
spacings (NNS), obtained from Si = Xi — and calculate the probability density function 



P(s) for a given sequence. It is useful to define the integrated probability function [20 
I(s) = Jq ds'P(s'). Notice that J(s) is nothing but the total number of NNS below a given 
s, and can be uniquely calculated, without any dependence on the specific binning used 
to calculate the typical histogram representations of P(s). In what follows, we use J(s) to 
make quantitative statements, but revert to showing the more conventional P{s). 

Figure 2 shows a sequence of N(E) staircase curves for different values of the background 
dielectric constant e, defined above. The sequences analyzed have all the same value of 
relative angular momentum, L Te \ = 3h, and total spin S = 3H/2, and for a magnetic field 
value of B = 3T. Increasing values of e would produce a progressively weaker value of the 
Coulomb interaction, and vice versa. We have varied e such that the interaction varies 
by up to a factor of 50. In the figure, the curve labeled 1 corresponds to the value of 
e = ei = 12.5, found in GaAs, where most quantum dots are defined. Increasing the 
effective Coulomb interaction (with e = ei/10), produces the smoothest staircase function 
(labeled 10), completely devoid of the harmonic oscillator 'steps' in weaker interactions (as 
seen here for relative interaction strength 1/5 and 1/3). One expects that for interaction 10, 
the level mixing produced would be substantial, resulting in a level structure well described 
by the RMT distributions, as seen in other systems. We find that this is not the case here. 

In fact, if one analyzes the appropriate histograms for the NNS distribution function P(s), 
the interaction strength shows its effect quite clearly, as shown in Fig. 3. There, we show 
the corresponding P(s) for three different interaction strength (or e) values. The anticipated 
Poisson distribution function one obtains for a 'generic' integrable system |||l| is similar to 
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the case for interaction 1/3. However, notice that this P here goes to zero even faster than 
the expected Poisson form, reflecting the non-generic character of the delta-function-like 
distribution of the pure harmonic oscillator system. Moreover, as the interaction becomes 
stronger, notice that the distribution functions have a maximum not at zero s, but rather 
at a finite value. This behavior is more in agreement with the RMT predictions, where P 
would be expected to be given by a GUE function (given the finite magnetic field), where 
P(s m 0) oc s 2 . One should notice, however, that in all cases we have studied, P(s) never 
fully reaches the anticipated GUE for a fully chaotic system. The reason for this lack of full 
crossover into the GUE distribution is perhaps associated with the mixed dynamics of the 
corresponding classical system. In such cases, it has been argued that the level distribution 
can be seen as a superposition of Poisson and GUE functions, with weights corresponding 
to the coverage in phase space for the integrable and non-integrable regions, respectively 
pOfl . At this point, however, we have not analyzed the classical system that directly models 
this quantum dot, and expect to report on this relation elsewhere. However, it may also 
be the case that this incomplete crossover to a GUE arises from the peculiar non-generic 
non-Poissonian distribution function of the integrable (non-interacting) system composed by 
overlapping harmonic oscillators, and/or remnant hidden symmetries. We have also fitted 
the /-functions to the well known Brody distribution, /erody = 1 — Pa exp(— s a+1 ), providing 
a (phenomenological) measure of the crossover (a = for Poisson, a = 1 for GOE). In Fig. 
2 and for unit interaction strength, we get a = 0.04, while a = 0.25 for interaction 10. 
Furthermore, the non- typical character of the level distribution in the harmonic oscilla- 



tor is also reflected in the magnetic field dependence pl| . If one studies P(s) for a given 



interaction strength but increasing magnetic fields, the system has a non-monotonic evolu- 



tion. Examples of this behavior are shown in Ref. [[H]. For low magnetic field (B = 0.3 T, 
and 'unit' interaction, with e = ex), P(s) is not the harmonic oscillator delta-like function, 
but is somewhat shifted towards a GUE (Brody-fit with a ~ 0.1). As the field increases, 
however, the P-distribution reaches a 'maximum' crossover (a ~ 0.2, for B = 1 T), before 
going back towards a more delta-like function, as seen for B = 10 T, which drops faster than 
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a Poissonian distribution. That this occurs is understandable, since for high magnetic fields 
one would expect to reach a regime where the interactions would be a weak perturbation 
(for a given value of e), and the spectrum would evolve towards a set of Landau levels. This 
non-monotonic behavior is seen in all the level sequences we have studied. 

IV. CONCLUSIONS 

We have shown that the interactions introduce classical non-integrability in a quantum- 
dot system, even if the geometry is integrable for the one-particle problem. Moreover, 
this behavior, given the Bohigas et al. conjecture 0, would be expected to be reflected 
in the accompanying level structure of the corresponding quantum mechanical version of 
the system. We have found, that as we introduce interactions, indeed the NNS fluctuation 
distribution function exhibits an apparent crossover towards one of the RMT functions, 
with a characteristic maximum away from zero spacing. This crossover, is found to be not 
complete, however, even for rather strong interaction, perhaps due to the mixed dynamics 
in the classical counterpart. The NNS distribution exhibits non-monotonic magnetic field 
dependence, associated with the fact that in that regime the particle interactions are but 
a weak perturbation of the Landau level spectrum. A number of theoretical questions 
still remain, however, as more direct comparison of the classical and quantum systems is 
needed, perhaps including a study of the associated wave functions and possible appearance 
of 'scars' ||. Nevertheless, the importance of interactions in the details of level structure 
and associated experiments, as discussed in the introduction, cannot be ignored. We hope 
that this, as well as semiclassical treatments (see Richter's article in these proceedings), 
motivate further understanding of this fascinating problem. 
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FIGURES 

FIG. 1. Poincare section of two-particles in a box with interaction given by A = 0.1. Total 
energy E = 1.5, and q\ = —1. 

FIG. 2. Typical staircase function for different interaction strengths. For interaction strength 
10, the harmonic oscillator structure is completely lost, while it is visible for weak interaction (1/3 
and 1/5). 

FIG. 3. NNS distribution function for varying interaction strength values (in units of the in- 
teraction in GaAs). Notice evolution towards a GUE distribution for stronger interaction. 
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